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Abstract 

Color superconductivity in QCD breaks the SU{2>) color gauge 
group down to SU (2), inducing masses in five of the eight gluons. This 
is a dynamical Higgs effect, in which the diquark condensate acts as 
the vacuum expectation value of a composite scalar field. In order to 
analyze this effect at low quark density, when gaps are large and gen- 
erated nonperturbatively, we use instanton-induced quark interactions 
augmented with gauge-invariant interactions between quarks and per- 
turbative gluons. The five gluon masses are found from the static 
limit of the relevant polarization operators, in which transversality 
is maintained via the Nambu-Goldstone modes of broken color sym- 
metry. Working in the microscopic theory we calculate these masses 
to one-loop order and estimate their density dependence. Finally, we 
speculate that the Meissner effect may postpone the onset of color su- 
perconductivity to higher matter density than estimated previously. 
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1 Introduction 



It has recently been demonstrated that the color symmetry of Quantum Chro- 
modynamics (QCD) might be spontaneously broken through the formation 
of a diquark condensate. This phase is metastable at zero and small matter 
density [|l| but becomes stable and replaces the usual chiral-broken phase at 
some critical density 0, |^. 

Condensation of diquarks is analogous to Cooper pairing of electrons in 
BCS theory, in which the massive photon is the microscopic embodiment of 
the Meissner effect. In QCD the gauge group is SU{3), but the consequences 
are similar: gauge bosons (gluons) become massive due to their interaction 
with the correlated fermions (quarks). Not every gluon species becomes 
massive, since the symmetry breaking pattern is SU{3) SU{2) and the 
three gluons of the residual SU{2) remain massless. Among the five massive 

modes, four are degenerate and the mass of the fifth is y^4/3 times that of the 
others. This is comparable to the symmetry-broken phase of the electroweak 
theory. Thus the direct analogy is to the Higgs mechanism, with the diquark 
field forming a composite scalar and its nonzero VEV generating a dynamical 
Higgs effect. 

In this paper this mechanism is analyzed microscopically. As the source 
of diquark condensation we will use the effective quark action obtained by 
averaging over instanton configurations |Q. In this approach to the QCD 
vacuum, all information from the large background gauge fields (instantons) 
is encoded in the non-local form of the 't Hooft interaction and the two 
instanton parameters, average size (p) and number density (N/V). This 
procedure replaces dynamical gluons with classical, nonperturbative field so- 
lutions, ignoring quantum gluonic fluctuations. The reintroduction of gluonic 
fluctuations will require a gauge-invariant perturbative modification of the 
't Hooft vertex. A direct consequence of gauge invariance is the transverse 
polarization operator, which we will compute explicitly. This is the main 
result of the paper as from it we obtain the gluon masses. 

Since the dynamical Higgs effect is interesting regardless of density, we 
will first consider the vacuum alone (chemical potential /i = 0), although in 
this case the superconductor is only metastable The gluon masses are 
found to be proportional to gA where A is the superconducting gap, the 
scale of which is a nonperturbative one set by the strength of the instan- 
ton background. Extending the formalism to finite chemical potential /i is 
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straightforward, given that the density dependence of the propagators and 
effective action are known 0]. However, in this paper we avoid exact cal- 
culations at /X 7^ 0. Instead, for simplicity we assume that the in-medium 
gluon masses are determined by the behavior of 5'A(/i) up to and somewhat 
above the critical value of /i ~ 300 MeV. At this density, the point at which 
the color superconductor becomes the stable phase, we find that the gluon 
masses are about 120 MeV. 

For gluons the color-breaking Meissner masses are a primary ramifica- 
tion of the superconducting state, however this is not the only effect of a 
quark medium. There will be additional contributions, among them a De- 
bye screening mass of the order of gfi, which do not break color symmetry. 
When the density becomes large, the Debye mass increases and instantons 
are screened out of the picture. Yet diquark condensation persists, now due 
to perturbative gluon exchange [Q, and at asymptotically large density the 
Meissner masses are also proportional to gfi So while one still has 

quark pairing, any matching between low and high density mechanisms re- 
mains unclear. 

In this paper we consider the case of two massless flavors, which is ex- 
pected to be relevant at finite-density chiral restoration given a relatively 
large strange quark mass p], ||. Instantons are Euclidean pseudoparticles, 
and thus all calculations will be in Euclidean space. 

In Section 2 we recall those features of the ordinary Higgs mechanism 
(based on elementary scalar fields) which will also be relevant for composite 
scalars. In Section 3 the instanton-induced action including perturbative 
gluons will be formulated. In Section 4 the diquark gap equation is reviewed, 
and in Section 5 the color current is determined. In Section 6 we sum a set 
of quark-quark correlation functions to recover the Nambu-Goldstone modes 
of the theory, a necessary exercise as these will mix with the longitudinal 
gluons. These ingredients are assembled into the gluon polarization operator 
in Section 7, which taken in the static limit corresponds to a mass. Section 8 
compares the result to chiral symmetry breaking, and in Section 9 this result 
is discussed and conclusions are drawn. 
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2 Higgs Mechanism 



The dynamical Higgs mechanism, though technically more involved, does 
not differ much from the ordinary one. Let us denote the quark bilinear 
combination 

0" = e^p,e^,e,, (^f ^ + (L ^ R)) , (1) 

where Greek letters denote color, f,g= 1, 2 flavor, and i,j = 1, 2 are spinor 
indices of the left (L) and right {R) components of the quark field ip. The 
resulting complex field 0" is a Lorentz scalar isoscalar field belonging to the 
fundamental representation of the 5^7(3) group. To support gauge invariance 
it must couple to the gauge potential via the covariant derivative, (V^)^ = 
dfj.S'^ — igA'^^{\°'/2)'p, in which the A*^ denote the eight Gell-Mann matrices. 
The kinetic energy term for the composite diquark field 0" is the usual 



(2) 



The only difference with the standard case of the elementary field is that 
there is, in principle, a common 'wave function renormalization' factor, Z^. 
For elementary fields Z^ = 1 aX the tree level, however, it deviates from unity 
even for the elementary field when one takes into account the perturbative 
virtual emission of particles. For the fields that are composite from the start 
there is no reason for to be unity. This quantity is a priori unknown and 
should be determined from a dynamical calculation, as will follow below. We 
stress, however, that the relative weights of the three terms in Eq. are 
fixed by gauge invariance. 

If the scalar field develops a nonzero VEV signaling the diquark con- 
densation, 

(r) = 2Ao5"^ (3) 

(it can be always arranged along the third color axis and made real), then 
gluons obtain a mass matrix 

a = 6 = 1,2,3 

Ml = 2g'Z^^/\l [X'X't = \ 2g^Z^'^l a = 6 = 4, 5, 6, 7 (4) 

a = 6 = 8. 
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Figure 1: Diagrams contributing to the gluon polarization operator in the 
Higgs effective theory. 



The symmetry brealcing pattern is, thus, SU{3) SU{2); three gluons 
corresponding to the unbroken SU{2) subgroup remain massless, four gluons 

obtain masses proportional the Higgs VEV, and the fifth gluon is 



'4/3 times 

as heavy. This relation will be, of course, reproduced for composite Higgs 
fields as well, as it follows from symmetry considerations alone. 

This simple elementary-Higgs model also provides an alternative way to 
find the gluon masses, which we will generalize to the case of the composite 
Higgs field. One can compute the gluon polarization operator H^^(g), as seen 
in Fig. |I[ In Fig. [^a we take the linear coupling of the gauge potential A^, 
as given by the second term in Eq. (^), and iterate it twice. The interme- 
diate state in this diagram is the would-be Nambu-Goldstone boson, which 
contributes to the polarization operator 



(5) 



This contribution is purely longitudinal. 

Fig. 03 is the contact term; it gives a Kronecker delta contribution: 



Combining the two we get a transverse polarization operator. 



(6) 



(7) 



The fact that it is transverse reflects the color current conservation, even in 
case of a broken symmetry. 
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In the case of a composite Higgs field the tree diagrams of Fig. |1] will be 
replaced by loop diagrams, formed by quarks, however, the general setting 
will be rather similar to this simple case: there will be a pole contribution 
from the would-be Goldstone bosons in the intermediate state, and there will 
be a 'contact' term. The resulting polarization operator will be transverse, 
provided one takes a conserved color current. 

3 Gauge Invariant Effective Action 

Following a long tradition of instanton phenomenology, we assume the fol- 
lowing: 

• Instantons are the dominant nonperturbative contribution to low-energy 
QCD, and 

• Low and high momentum scales are safely separable for quark and 
gluon fields. 

These assumptions are supported a posteriorihj instanton-based phenomenol- 
ogy of the vacuum, as reviewed in Ref . []TU| . Instantons break chiral symmetry 



spontaneously and axial U{1) anomalously, both at a satisfactory magnitude. 
Instanton-induced interactions also provide the necessary qq attraction for a 
diquark condensate, nonperturbatively. 

Here we use a formulation which not only relies on the separation of quark 
zero and free modes, but also explicitly includes both instantons and pertur- 
bative gluons. Gluons are separated into classical instantons and quantum 
corrections, which we write as 

= E4'H^) + E<na^) + «.(:^), (8) 

where the sums are over all instantons and anti-instantons, each given by 
the 't Hooft solution in the singular gauge [^. Although there is a certain 
distribution of instanton sizes (see, for example, [jl2|) we simply use the 
average value, p ~ 1/3 fm, in all calculations. 

The main idea of this paper is to examine the effects of color superconduc- 
tivity on the gluonic excitations above the instanton vacuum, which means 
we must retain and analyze the a^. Hereafter, the term 'gluon' will refer to 
the perturbative gauge fluctuation above the instanton background, a'^. 
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The origin the 't Hooft interaction is explained in the hterature |TT[ and 
will not be repeated here. To summarize, the low momentum quarks are 
approximated by zero mode solutions in the presence of one instanton. Av- 
eraging over the instanton ensemble generates a vertex for dynamical quarks, 
mediated by the zero modes, which can be treated perturbatively when the 
instanton liquid is reasonably dilute. The relevant small parameter is the 
ratio of average instanton size (p) to the average inter-instanton spacing (R). 



Prom phenomenological [|T3| , variational [|T4[ , and lattice calculations one 
finds 

14- 

More details follow in Section 4. 

We now consider the effective action itself. In the quark sector it is a non- 
local one of four-fermion operators, since quarks are connected to instantons 
via the quark zero modes whose spatial extent is of the order of the instanton 
size, p. Specifically, one has 0]: 

Nf 

SiNT = X jdUd^z md^Xf dSf ^jixfW^Xf - z,U)^{yf - z,UWijf{yf)l 

(10) 

where U is the instanton's 2 x Nc color/spin orientation matrix and z is its 
position. The ^{x) is the zero mode solution for fermions in the field of one 
instanton; in general it depends on the chemical potential. At /i = its 
Fourier transform is the form factor, 



f{p) = 2x 



Ii{x)Ko{x) - h{x)Ki{x) + -h{x)Ki{x) 

X 



(11) 

x=pp/2 



such that /(O) = 1. At /i 7^ the form factor is also known explicitly 0. 

Because of non- locality Eq. (|I0|) is not gauge invariant. When calculating 
bulk vacuum properties one needn't worry about quantum corrections as 
observables are seldom sensitive to them, however in this case the dependence 
is crucial. A limited literature does exist in which non-local interactions 
are modified to be gauge invariant |]15|, [1^, |18| and we follow the same 
procedure here. In particular, we will minimally modify the effective action 
( p!OD to suit the present needs. 

We will take the non-local four-fermion interaction as the starting point. 
It then becomes a matter of multiplying each quark operator by a path- 
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ordered exponential in the background of the perturbative gluon field a^, 
replacing as: 



ip^x) — > ip{x)W{x, z), 
W{x,y) = Ve^vii^-jjs.alX^y (12) 

where is the perturbative field. As has been pointed out in the cited 
works, the choice of path integrated over is not unique. Yet as long as these 
factors transform as 

W{x,y)^U{x)W{x,y)U\y) (13) 

by virtue of path ordering, the action will be gauge invariant. This remains 
true when, as in this case, the interaction involves an explicit color average 
since overall color is conserved. We are also concerned only with the static 
limit, — > 0, in which results are independent of any particular choice of 
path. 

Thus we write the modified interaction as a product over fiavors, 



N 



SiNT = / II n^-^/ ^"^yf 

x^]{xf)W{xf, z)^^{xf - z)^\yf - z)^W{yf, z)^ljf{yf)\ , (14) 

noting that this includes terms of all orders in ga^^. Calculating the gluon 
polarization operator will require the linear and quadratic contributions. 

To this interaction term one must add the usual quark kinetic term, min- 
imally modified to preserve gauge invariance: 

Skin = J d'x ^h,{td, + |A"a;)^. (15) 

Consequently, the color current obtained from the variation of the action in 
respect to a° will have two contributions: one is the ordinary one arising 
from the minimal coupling (^) and the other arising from the non-local 
interaction term (0). 
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4 Diquark Condensation 



In this section some details of the color superconductor phase are reviewed. 
All expressions are Euclidean and the notation generally follows that of 
Ref. the reader is referred to this reference for details specific to this 
particular approach. The spontaneous breaking of chiral symmetry was also 
a central concern in that work, whereas here we are considering restored chi- 
ral symmetry. The main result of Refs. [IS, 0, 21, 22] is a competition 



between chiral and diquark condensates at zero temperature and nonzero 
quark chemical potential. In all cases, there is a phase transition from a low- 
density phase of spontaneously broken chiral symmetry to a high-density one 
of color superconductivity. 

While technically involved, our previous results arise from an effective 
four-quark interaction which allows for pairing of quarks as in the BCS the- 
ory. The instanton model retains more of QCD's features than ad hoc models, 
notably the anomalous breaking of axial U{1) and transmutation of dimen- 
sions, and is constructed in a more systematic way from that underlying 
theory. 

A prominent feature of the instanton approach used here is the deter- 
mination of the effective four-fermion coupling constant A. This constant is 
determined by a saddle-point evaluation and proves to be nonlinearly depen- 
dent on the background instanton density, N/V . A gap equation obtained 
from a set of Schwinger- Dyson- Gorkov equations are solved to first order in 
A, done self-consistently to determine the quark pairing gap, Aq. Standard 
quark propagators are necessarily split by a nonzero Aq, since this introduces 
a color bias and quarks are no longer color degenerate. This detail will be 
crucial when quark loops are computed in Section 7. 

Since this paper deals only with a phase of chiral- symmetric, diquark 
condensation the Schwinger- Dyson- Gorkov equations and corresponding di- 
agrams are the simple ones of Fig. 0. We define the quark propagators as 



I'^f'^HvU^ - <! ^l^P^^P'^^ a, 13 = 1,2 

[P)%pj[P)) - i SfS'^S2{p)i a,P = 3 ' ^^^^ 



and the anomalous Gorkov propagator as 
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Figure 2: Schwinger- Dyson- Gorkov diagrams for normal and anomalous 
quark propagators. The labels refer to the functions in the text. 



In these expressions indices / and g refer to flavor, i and j to spin, the Greek 
letters to color, and ipL,R are chiral spinors. Written in the chiral L, R basis, 
the 4x4 propagator Si{p) is of the form: 



Slip) 



z{p)So{py 

Z{p)So{p)- 



(18) 



while 5*2 (p) is the free propagator for color 3 and hence identical to the above 
with the function Z{p) absent. The notation is = x^a^, where the 2x2 
matrices cr^ = (±i(T, 1) decompose the Dirac matrices into chiral components. 
The off-diagonal, bare propagator is therefore written So(p)^ = 

With these definitions, we have the scalar Schwinger-Dyson-Gorkov equa- 
tions 

Z{p) = l-A{p)F{p) 

A(p)Z(p) , , 

F{P) = \ • 19 

The momentum dependence of the gap, 

A(p) = Ao/b)^ (20) 

is given by the instanton- induced form factor, Eq. (|lT]), which suppresses the 
interaction beyond p > l/p ~ 600 MeV. This pair of equations leads to a 
gap equation for Aq, 

- ^77^/7^^A(p)F(p) 



N,{N, - 1) J (27r) 
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2A r d^p Ajpy 
N,{N, -I) J (27r)4 p2 + A(p)2 ' ^ > 

which must be self-consistently solved with A. This couphng constant is in 
turn determined by a saddle-point integration (exactly in the thermodynamic 
limit N,V ^ oo), 

A - ^/^ ^0 • (22) 
Combining these two equations A may be eliminated and we obtain 

^^\V) 7 (27r)4 p2 + A(p)2 • l^^J 

The scale of the gap, therefore, is set by the instanton density N/V. To this 
order N/V remains at its vacuum value, and although it will be affected by 
the finite density of quarks this is an 0{X) correction to the instanton weight 
and not considered here. 

Solved numerically, the gap is Aq — 400 MeV in vacuum and drops to 
Aq — 200 MeV at the phase transition from chiral broken to color supercon- 
ducting matter, as detailed in Ref. 



5 Color Current 

A color-conserving Noether current naturally follows from the modified inter- 
action, including contributions from both Skin 0T5| ) and Sjnt (0)- Along 
with the standard quark-gluon coupling piece. 



d^p 
(2^ 



i'R{p)^(^/.i'L{p + q) + i'L{p)^(^t^R{p + q) 



(24) 



the current now includes a four-quark coupling to the gluon as shown in 
Fig. and written in terms of four-momentum: 



5S 



INT 



aj=0 



In Ref. PI Aq was defined as half this paper's (and the more standard) definition. 
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Figure 3: Quark-gluon vertices of order g (left) and gX (right) interactions. 

_ ^9192 J j^^-ixf{pf-p'f)+iyf{kf+k'f)-iz-(j>'f+k'f) 



X 



i , 1 1 

ds,e-'^-'^ljif^{p,)M{p[, k[)X'rL\h)ijif^{p,)M{p'^, KWL\k2) 
+{L^R). (25) 
Indices on quarks denote chirality and flavor. The measure is 

T. ^rM4 ^TT^^^ ,4 d^Pfd*p'fd*kfd% 

V^dUd z [[ d Xfd Vf j^^^ (26) 

and the form factors lie in the color /spin matrices 



A^(P, m = U^e%iUl^f{p)f{k). (27) 

The U are 2x A^c color orientation matrices, averaged in each vertex. Through 
use of the Dirac equation for ip , one can explicitly verify that q^J^{q) — for 
J^l = jfi + Jfi- While this condition would remain satisfied with a transverse 
addition to the color current, no such addition is motivated here. 

In practice we are interested in pairing off two of the four quark legs of 
the second vertex. Since we are considering a phase where chiral symmetry 
is unbroken, a chirality-violating {ij^ip) loop cannot contribute. Thus we pair 
either (V'V') or {ip'^ip^) and obtain the effective current, 

^0 f..d%d%d'p' 



(27r) 
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X 



^-ix.{p,-p')-iz-(j>'+p2) r dSf,e-''^-'ij[{pi)X''X{p2,p')^i{p2) 

Jx 
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+ {L ^ R) . (28) 

where here the flavor /color /spin structure is 

U{p, p')//''""''"' = ^e^^'e""'3e"7(p)/(p') • (29) 

The two terms in Eq. ( pH]) correspond to pairing either two incoming or two 
outgoing quark legs of Fig. ^ 

Gluon mass terms are found in the — > limit. When the gluon couples 
to nonsingular composite quark modes one can set strictly to zero and the 
path integrals simplify to 



= e-'P^(^-"^)(z-x,)^ 

= -ie-'P''<'-'''^^. (30) 
dp'if. 

This substitution, the result of which clearly depends only on the endpoints 
of the path rather than any particular choice of path, leads to differen- 
tiation of the form factor. On the other hand, the vertex can also couple 
gluons to Nambu-Goldstone modes which are singular as 0. Their 
behavior must be countered by expanding J^(g) to order q^, but before this 
complication arises the Nambu-Goldstone modes must be specified. 



6 Nambu-Goldstone Modes 

A symmetry has been spontaneously broken and Nambu-Goldstone modes 
are certain to follow. With SU{3) being broken to SU{2) there are five 
massless modes and, directly analogous to the Higgs mechanism, they do not 
become additional degrees of freedom. Instead they mix with and are incor- 
porated into the five massive gluons, thereby relevant in the color Meissner 
effect. 

Since the associated condensate is (qq), these massless modes will be 
quark bilinears. To determine their quantum numbers one need only perform 
a gauge rotation on the diquark condensate and catalog the five orthogonal 
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Figure 4: Bethe-Salpeter-Gorkov diagrams for summation of the quark-quark 
correlation functions. 

correlators which appear. Recalling the condensate direction as chosen in 
Eq.®, 

Z6.,6;,(^fV^f)~Ao, (31) 

we can list the five diquark operators which couple to the Nambu-Goldstone 
excitations: 

^^fg^kl^l'^^^^L^ ^fg^kli^L^^^'^'^L ■ (32) 

Propagators for these quark-quark modes will exhibit a simple pole at = 0, 
behavior recovered by computing the corresponding correlation functions. In 
order to obtain this pole we must sum s and u channel contributions to all or- 
ders in the four-quark coupling A 0. Since connecting these modes to the gluon 
propagator will be itself of order g"^, we do not include any gluonic corrections 
to the instanton vertex. With both standard and anomalous quark propa- 
gators at our disposal we obtain the set of coupled Bethe-Salpeter-Gorkov 
equations diagramed in Fig. ^. The diagram shows not only standard two- 
body propagators, denoted F, but also its anomalous analog, Q. It is easy 
to see that Q will vanish when any of the external lines are quarks of color 
3, since the vertices conserve color and the internal (Gorkov) propagators 
involve only colors 1 and 2. 

Knowing the quantum numbers of the Nambu-Goldstone modes, we can 
immediately write down ansatze for the four-point functions which will be 

related summation could be done for t channel diquark exchange as well, however 
these do not contribute the dominant l/q"^ behavior in the limit of zero transfer momentum. 
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required when we later compute the gluon polarization operator: 



N]{N, - 1) 



X (27r)^5^ ip + p'-k- k') fip)fip')fik)fik')Tf^,ip + p') (33) 



1/3/3' 



X (27,^5' ip + p' + k + k') fip)fip')fik)fik')niUp + p') . (34) 



Here, [ea]^^ = eso,/?, and the x refer to L or i? chirality with the chiral 
substructures defined as 



^RL ^RR 



(35) 



in which color indices have been suppressed. It is clear from their definitions 
that Tlr = Trl, Tll = ^rr, ^lr = ^rl, and Qll = ^rr- 

Written in terms of these chiral elements, the two diagrams of Fig. ^ 
correspond to the following set of equations: 





= A 






= A 








= A 


'ciT'tMUq) + c''A{q)Uq) 








= A 


'c^rf^{q)Uq) + c^,Qt,{q)Uq) 


5 





(36) 

where no sums are implied on color indices a, b. The adjoint color vectors 
{ci} determine the internal propagators and are found to be 



Cl = (0,0,0,1,1,1,1,0) 
C2 = (1,1,1,0,0,0,0,1) 
C3 = (-1,-1,-1,0,0,0,0,1) 

The integrals which result from these loops are: 



(37) 



d'^p p+ ■ p- 



J (27r)4 pIpI 



.[A(p+)A(p_)] [Z(p+) + Z(p_)] 
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where we have defined (j9±)/x = P^J,^\(l^l and 

The set of Bethe-Salpeter-Gorkov equations (^) can easily be solved for 
correlations functions of all adjoint colors. For colors a = b = 4 ... 8, in 



direct analogy to similar calculations of pions in the instanton vacuum 
the = pole arises due to a cancellation in the denominator which is a 
direct consequence of the gap equation. The calculation follows identically 
for a = 4,5, 6, and 7, where we find 



After writing the gap equation ( psf ) in these terms, 

l = AJi(0), (41) 

we have 



where 



ZT^ = -2 



dq 



g2=0 

rfV rA(p)2-iA(p)pA'(p) + ip2A'(p)2 



AIJ (27r)4[ p2[p2 + A(p)2] 

iA(p)2[A(p)-pA'(p)] 



p2|'p2 _j_ A(p) 



212 



(43) 



The primes denote differentiation with respect to p. The value of Z^ will 
change with density, since not only does it depend on Aq but the integrand 
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will exhibit functional dependence on /x. In vacuum numerical evaluation 
gives Z^^ = 8.07 x 10"^ 

For the F^^ and solving the coupled equations becomes more involved, 
since all diagrams in Fig. H are present. One finds 



A 



l-A2[X2(g2)+X3(g2)]2 

= (44) 

where the second line is obtain by making use of the Schwinger-Dyson-Gorkov 
equations, Eqs. (|19]), in the limit of small g^. This completes the set of 
Nambu-Goldstone modes. One can naturally compute the correlation func- 
tion for the additional diquark correlators with a = 6 = 1, 2, 3, and find a 
crucial sign difference in the combination of integrals: 



A 



1 



g2 + m{qy 



(45) 



The mass can be determined from an 0{q^) expansion of the denominator, 
however this will not be done here as for our purposes it is enough to verify 
the a = 1, 2, and 3 modes are indeed massive. 



7 Gluon Polarization Operator 

With the gapped quark propagators, conserved current interactions, and 
composite modes defined in the previous sections, we now compute the lead- 
ing modification to the gluon polarization operator. All diagrams prove to 
be color diagonal, and so we write 

nJ(g')=5"U.(g'). (46) 

Our interest is the static limit, q^ — * 0, which may be considered an effective 
mass. 

In the presence of a color-3, scalar diquark the gluons are divided into 
three classes. Gluons of adjoint colors 1, 2, and 3 belong to the residual SU{2) 
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(5a) (5b) (5c) 






(5d) (5e) (5f) 

Figure 5: Quark loop diagrams contributing to the gluon polarization oper- 
ator. 

gauge group and as such remain massless. Gluons 4, 5, 6, and 7 couple one 
gapped quark (of fundamental color 1 or 2) with the ungapped species and 
share a degenerate mass. Gluon 8, diagonal in fundamental color, obtains a 
mass proportional to the previous four. One polarization operator from each 
class will be explained here to avoid unnecessary repetition. All possible 
contributions to order are diagramed in Figs. |^ and ^ and, depending on 
the gluonic species, some of these diagrams vanish and others combine to 
cancel in the static limit. 

We begin with the case of gluons 4-7, considering corrections to (a^a^). 
Diagrams (^), (^), and those involving (^) require pairs of Gorkov prop- 
agators and thus vanish since the quarks of color 3 (to which these gluons 
couple) cannot propagate anomalously. There are four graphs which consti- 
tute the 'contact' term proportional to 6^^- 

Diagram (^) is a standard loop, where one of the quarks is gapped (color 
1 in the case of gluon 4) and the other not (color 3). After subtracting off 
the vacuum part of this diagram, which remains a concern of gluon renor- 
malization and is not relevant to the Meissner mass, we find 

n<?>tf^0)^ViVA./(g,^55^^^. (47) 

The remaining integral is finite due to an power-law cut-off in the the function 
A(p) arising from the finite size of instantons. 

Additional diagrams arise from the modified Sjnt- Not only does this 
generate additional current interactions (pOf), but the interaction ( p!^ itself 
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contains a contact contribution to the gluon two-point function. This, di- 
agram (||c), is the second variation of the action (|l^) with respect to the 
fourth gluon field, 



6'S 



INT 



(48) 

g2=0 



^ 6af,iq)6atiq) 
Evaluated to order g"^ it is 

Diagrams (||d) and (|e), constructed with the additional current piece, are 

Eqs. (|47|), (^), and (0) comprise the microscopic equivalent of the Higgs 
contact term (^), and their sum is 

n?j(g^-o) = 



^^^^'^ I (2^ P'[p' + A{py] • ^^^^ 

The integral can be trivially rewritten and then manipulated to yield 
rfV rA(p)2-iA(p)pA'(p) 



(27r)4 [ p^[p^ + A{p) 
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A(p)pA'(p) A(p)p2A"(p) 



4p2 [p2 + A (p) 2] 4p2 [p2 + A (p) ^ 



Z7^A^ . (52) 



The final equality is achieved by integrating the second and third terms by 
parts. 

The Nambu-Goldstone modes couple to the gluons as in Fig. ^ The con- 
struction [(|^a) + (PD)]r^^(g^) [(|a) + (|b)] supplies the q^qu/q^ piece to ensure 
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(6a) 



(6b) 



(6c) 



Figure 6: Diagrams coupling the Nambu-Goldstone modes to gluons. Their 
sum squared contributes to the polarization operator. 



transversality: 



0) 



(27r) 



p2 Jp2 -|- 



A(p)2 - ^A{p)pA'{p) + \p'^A'{pf 



1 2 



p2 [p2 + A(p)2 



q2 



(53) 



We have now accounted for all contributions to gluons 4, 5, 6, and 7. 

Analysis of the eighth gluon follows in a similar fashion, although with 
terms from every diagram. A superficial difference lies in the factors arising 
from the elements which lead to a polarization 4/3 times the previous 
result. More subtle is the combination of diagrams (^) and (^). For this 



sums to 4/3 times Eqs. (^7[) , whereas for a^'^'^ they cancel one another. This 
cancellation is only manifest to 0{A'^ here, due to the fact that the diagrams 
are constructed with a pair of one-loop, resummed quark propagators as 
determined in Section 3. This generates a A{pY term in each integrand 
denominator, in essence including higher-order terms in the perturbative 
expansion (in A) which violate gauge invariance. After the vacuum pieces 
are subtracted, these diagrams contribute the following to a^: 



n(5;^|)(,^^0) = Ig'NfA^^j 



p^ 



d^p r 

(2^\F+AM 
2A{pf 1 1 
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[p2 + A{p) 



212 
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4 2.. ^2. f d'p \ A{pf 



(54) 

For gluons 1, 2, and 3, the sign of the second term is changed and the 
corresponding integral is of order Aq, and thus this quantity vanishes to the 
order to which we know quark propagators, Aq . The same follows for pairs 
of diagrams with similar construction, such as (§e) and (^). 

By combining Eqs. (|46|), (pT]), ([5^), and ( |53| ) we arrive at the satisfyingly 
compact expressions 

r a, 6 = 1,2,3 

Uf,{q' 0) = (5,, - ^) a, 6 = 4, 5, 6, 7 (55) 

[ -Ig^NfAlZ^'S^'lS^.-^-f^) a,b = 8. 

Transversality requires that the contact term be proportional to the wave 
function renormalization of the Nambu-Goldstone modes and this result, a 
Ward Identity for color superconductivity, is recovered here. 

Finally, to determine a numerical value for the masses we must fix the 
coupling constant g. Evaluating in the instanton vacuum, one finds the large 



finite action 



^0 = ^^12, (56) 

g 

or (yf ~ 2.6 and the perturbative expansion parameter = g^/4:iT = 0.54. 
The gluon masses squared are thus 

a = 1,2,3 

= { 2g^AlZ^^ ~ (150 MeV)^ a = 4, 5, 6, 7 (57) 

1^ 



]g^AlZ>^ ^ {UbMeVy a 



These masses apply to the vacuum, /i = 0. 

In order to estimate the finite-density behavior of the Meissner mass, we 
can simply take the values of Aq for a given /i from the results of Ref. 
As detailed in that paper, the instanton form factor (0) becomes density 
dependent and thus A(p) should be replaced by a comphcated A(p4, |p|,/i). 
However, the changes in Z^ arising from the finite- modifications of the form 
factors are minor compared compared to the changes in the gap magnitude. 
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100 200 300 400 500 600 
|x (MeV) 



Figure 7: Scaling of M, , and Aq with quark chemical potential /i. Each 
quantity is shown relative to its vacuum value. For reference, diquark con- 
densation is only stable in quark matter at and above /i ^ 300 MeV. 



Aq; we numerically estimate this correction to be about 3%. For simplicity we 
therefore considered only the scaling from Ao(/x) (taken from previous work 
[§]). Each non-zero gluon mass is proportional to AqZ^^''^ (see Eq. (pTf)) 
and therefore all will scale identically with density. The resulting gluon 
mass M(/i), as well as that of the renormalization constant Z0(/i)~^/^, is 
shown in units of its vacuum value in Fig. |^. Note that although rises 
with increasing quark density, this effect is not sufficient to overcome the 
falling gap and the masses continuously decrease. The first point of physical 
relevance would occur around /i ~ 300 MeV, the common prediction for 
chiral restoration to a color superconductor. Results for matter at lower 
densities correspond to an unstable solution and are only of academic interest 

a ra, ^,i2Ti,Ei. 
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8 Similarity with Chiral Symmetry Breaking 



Apart from the gauge coupling g the gluon masses (|57D are determined by the 
combination = 2AqZ^^. This quantity is the analog of the constant 
in the chiral-broken phase, both in its physical meaning and algebraically. 

In the chiral-broken phase the Nambu-Goldstone bosons are pions. The 
correlation function of the axial current, j^g = '?/'^7^75r^'?/', in the massless 
quark limit has the transverse form 

= n^.^(g) = -f^J""^ - ^) . (58) 

The transversality of 11^,^ is the consequence of the conservation of the axial 
current; the pole is due to the pion in the intermediate state. Were there 
gauge bosons coupled to the quark axial current their mass would be equal 
to Ft, multiplied by the corresponding gauge coupling. In case of diquark 
(vs. quark-antiquark) condensation the relevant currents are color ones, and 
we obtain a similar form for the correlation function of two color currents, 
Eq. (^3]). One needs only to multiply Fgg by the gauge coupling to deduce 
the gluon mass. 

If instantons are dilute, the leading contribution to the Kronecker part of 
the polarization operator arises from the Ti.^^^^ piece, Eq. (^Tj). It is the only 
contribution which diverges logarithmically if one neglects the momentum 
dependence of the gap A(p). For the same reason the leading contribution 
to the Z^^ and F^^ in the dilute limit comes from the pieces not containing 
the derivatives A'(p). One has therefore in the dilute limit: 

p2 ^ 2 / - — ~ ^ loff f5Ql 

J (27r)4 p2[p2 + ^(p)2] - 8^2 p2 

where Aq is the superconducting gap at zero momentum. It follows from the 
gap equation (|23|) that Aq ~ np^N/V = np/R"^. Similarly, the axial corre- 
lation function (|58| ) computed in Ref. |^ gives (in the same approximation) 



~ 4N fj^^dlit_ ^ 1^ log ^ (60) 

where M(p) is the dynamical quark mass (the chiral gap) whose value at zero 
momentum is determined from a corresponding gap equation to be Mq ~ 

t^p^nJvn,. 
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We see, thus, that not only have the Meissner mass and the F.^ constant 
analogous meaning, but their algebraic structure is quite similar. One ex- 
pects, therefore, that the numerical value for the Meissner mass is of the order 
of F^, and this expectation is confirmed by an exact numerical calculation of 
the previous section. 



We have analyzed the problem of spontaneous gauge symmetry breaking 
brought about by a diquark condensate. Since the broken symmetry is con- 
tinuous and gauged the resulting Nambu-Goldstone modes do not remain in 
the spectra, instead mixing with the longitudinal components of the gauge 
fields to produce massive gauge bosons. This, a dynamical Higgs mechanism, 
can be called the color Meissner effect in the context of color superconduc- 
tivity. 

To reveal the gauge boson masses mathematically one has to compute 
the polarization operator, which ought to be transverse. 



where the massless pole l/q"^ is the manifestation of the Nambu-Goldstone 
intermediate state. The coefficient, gives the mass of the gauge boson. 
In this paper, we exphcitly solved this problem for the case of diquark con- 
densation as induced by the instanton background, the effective action from 
which we found necessary to modify in order to maintain a conserved color 
current. 

Through computing the gluon polarization operators to order g'^ we find 
the effective gluon masses to be on the order of the diquark gap. The three 

gluons comprising the residual SU{2) group remain massless and hence a 
quark-gluon medium would become color-biased in such a phase. The anal- 
ysis here was done in the limit of zero temperature and, initially, vanishing 
chemical potential //., though strictly speaking at zero n the color supercon- 
ductor is only metastable with the ground state being the usual chiral-broken 
phase. We then estimated finite-density dependence of the calculated quan- 
tities. At the critical density, where one expects the phase transition to the 
color superconducting phase, we deduce that the Meissner masses of the 4*'*, 



9 Discussion and Conclusions 




(61) 
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5^^, 6^^, and 7*^ gluon are about 120 MeV and the 8^^ gluon has a mass 
about 140 MeV. These quantities would be of physical relevance should a 
low temperature, high density region become experimentally accessible. At 
a chemical potential low enough to leave the instanton background approxi- 
mately unchanged (/x < 0.6 GeV), the instanton effects analyzed here would 
still be present and likely dominant. 

In computing the coupling of the Nambu-Goldstone modes to gluons we 
have established that their mixing is described by an effective theory in which 
the composite diquark is replaced by a complex scalar in the fundamental 
representation of the gauge group, 0". The effective Lagrangian coincides 
with that of elementary scalar field covariantly coupled to gluons, aside from 
an overall factor which is the 'wave function renormalization' of the com- 
posite scalar field. There are no a priori reasons for this factor to be close to 
unity (as it is in the case of a weakly-coupled elementary Higgs field), and 
indeed we find a substantial deviation from unity. Meanwhile, it is a crucial 
factor for the estimate of the Meissner mass. 

If electromagnetic interactions were taken into account, the gluon will 
invariably mix with the (massless) photon. This mixing, estimated from 
general arguments to be rather small [^], reorganizes the fields into a mas- 
sive 'new gluon' and massless 'new photon'. Given that the gluonic sector 
alone can be recast as a Yang-Mills-Higgs theory, this additional result would 
complete the analogy to the electroweak sector of the Standard Model. 

Finally, we would like to comment on the phase transition between ordi- 
nary chiral-broken phase and the color superconductivity. The point is that 
all estimates existing in the literature (see Ref. and references therein) 
indicate that this phase transition happens alarmingly 'early': taken liter- 
ary, the claim is that the interior region of a heavy nucleus is actually in a 
'boiling' state. However, those estimates generally neglect the influence of 
the dense medium on the gluonic background fluctuations which induce the 
color superconductivity itself. The Meissner mass of about 150 MeV found 
here is a large quantity and, together with the Debye mass and other effects, 
will suppress instantons. Therefore, by taking into account back-influence 
effects one has a chance to 'save' ordinary nuclear matter from a premature 
phase transition by moving that transition to higher densities. 

Although the instanton density is not expected to change significantly for 
any chemical potential below the inverse instanton size (600 MeV), at and 
beyond this point perturbative Meissner and screening masses will become 
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increasingly important. At asymptotically large densities perturbative gluon 
exchange is the source of diquark condensation 0, ^, ^ ^ whereas 
in this work we have been concerned with the low density, nonperturbative 
regime. Determining the behavior at a moderate density scale - which would 
be an interpolation between the two - is necessary if one wishes to confidently 
consider signals of color superconductivity in any potentially realizable situ- 
ation, be it in some future heavy-ion collider or a neutron star. 
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